Non-Hermitian Spectra and Anderson Localization 



00 

o 
o 

> 

o 
O 

m 



43 : 

Oh- 

> . 



> 
od 

o 

oo 
O 



Luca G. Molinari 
Dipartimento di Fisica and I. N.F.N, sezione di Milano 
Via Celoria 16, 20133 Milano, Italy 
E-mail: luca. molinariQrni. infn. it 
(Dated: October 2008) 

The spectrum of exponents of the transfer matrix provides the localization lengths of Anderson's 
model for a particle in a lattice with disordered potential. I show that a duality identity for deter- 
minants and Jensen's identity for subharmonic functions, give a formula for the spectrum in terms 
of eigenvalues of the Hamiltonian with non-Hermitian boundary conditions. The formula is exact; 
it involves an average over a Bloch phase, rather than disorder. A preliminary investigation of 
non-Hermitian spectra of Anderson's model in D=l,2 and on the smallest exponent is presented. 

PACS numbers: 71. 23. An (theories and models, localized states), 02. 20. -a (Matrix theory) 



I. INTRODUCTION 

Several models in physics are described by matrices 
with banded or block-tridiagonal structure. Examples 
are the Laplacian matrix, the Anderson Hamiltonian for 
transport in a lattice with random impurities, band ran- 
dom matrices, tight binding models in condensed matter 
and chemistry. The matrix structure reproduces that of 
a system consisting of a chain of units with same num- 
ber of internal states, with nearest neighbors interaction. 
Finite size effects are often dealt with by imposing peri- 
odicity; the limit of large number of units is eventually 
taken. 

The matrix structure calls for a transfer matrix de- 
scription of the eigenstates, and the spectrum of expo- 
nents of the transfer matrix describes the decay lengths 
of the eigenstates. For the Anderson model or band ran- 
dom matrices, most of the knowledge on Lyapunov spec- 
tra relies on numerical computations. 

With great generality, I showed that an analytic tool 
to access the decay lengths is a duality relation, that 
connects the spectrum of a non-Hermitian extension of 
the block-tridiagonal matrix, with the spectrum of the 
related transfer matrix2u2i2is. The extension arises by 
mere generalization of boundary conditions for the eigen- 
states. If Uk specifies the state of a unit of the chain 
(k = l,...,n), the boundary conditions (b.c.) are 
parametrized by a complex number z: 



U n+ i = Z Ui, Uq = —u r . 

z n 



(1) 



This destroys Hermiticity of the Hamiltonian matrix, but 
enlights a nice property of the transfer matrix: 



T(e 



Ul 




Ul 




= z n 











(2) 



The ensuing spectral duality and Jensen's identity for 
subharmonic functions allow to evaluate the counting 
function of exponents. This paper is intended to intro- 
duce the theory and explore its application to the long- 
studied problem of Anderson's localization. In section 



II the duality relation is reviewed and the main formula 
(|12p for the exponents is obtained from Jensen's theorem. 
The theory can be extended to include the spectrum of 
the (Lyapunov) exponents of the matrix T^T, by con- 
structing a corresponding non-Hermitian block tridiago- 
nal matrix, twice the size of the original Hamiltonian ma- 
trix. In section III a preliminary study of the eigenvalues 
of non Hermitian Hamiltonian matrices in D=l and D—2 
is made, with the purpose of illustrating the duality. The 
spectral formula is used to evaluate the smallest exponent 
imin, in a regime where the eigenvalues are already all 
complex. 

Because of their relevance in mathematics, numerical 
analysis and physics, block tridiagonal matrices are an 
active area of research^! 7 -. This work extends in a new 
perspective the work by Hatano and Nelson^ which, to- 
gether with the works by Feinberg and ZeeS, started an 
interest for non-Hermitian matrix models in physics. 



II. THEORY 

A. Transfer Matrix 

Consider the following block tridiagonal matrix with 
corners, of size nm x nm, 



H 



A x B 1 
B\ ■. 

B n 



4-i 



An 



(3) 



The blocks have size m x m: Bk are complex matrices 
with det Bk ^ 0, Ak are Hermitian matrices. To the 
matrix el nm — H there corresponds the transfer matrix 1 



T(e) 



B-\eI m -A n ) -B^B\_ X 

Irn 



(4) 
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I m is the to x to identity matrix. The transfer matrix is 
so named because it transforms the eigenvalue equation 
Hu = eu into a relation for the end-components of the 
vector u — (ui, . . . , 



T(e) 







U n +1 









(5) 



The corners imply a condition of periodicity uq — u n and 
Un+i = that can be used to obtain the eigenvalue 
e in alternative to diagonalization of H. By comparing 
eqs.Q and ([5]) one arrives at the main point: to study the 
spectrum of T(e), one must impose the generalized b.c. 

We thus introduce an instrumental non-Hermitian 
matrix depending on a parameter z (0 < argz < -) 



H(z n ) = 



At Bi 
B\ ■■• 

z n B„ 



B 



-Bi 



B n -i 



(6) 



The matrix is Hermitian for Bloch b.c. (\z\ = 1) but, for 
the purpose of studying the spectrum of T(e), it will be 
considered for z 6 Cq. The matrix can be brought by 
similarity to the balanced form H^z) = Z^ 1 H(z n )Z , 



H b (z) = 



A x zB x 
\B\ 



zB r , 



L -B^ 



zB n ^ 
iRt A 

z n n-l A n 



(7) 



by means of the block diagonal matrix Z with blocks 
{zl m , . . . , z n I m }. Therefore, no site of the chain is priv- 
ileged. While the matrix Hb(z) does change if arg z is 
increased by 2tt /n, its eigenvalues do not. 



B. Symplectic properties and exponents 

The following relations hold for the transfer matrix 
T(e*)tS n T(e) = E„, S„ = 
T(e)X- l T(e*)1 = E" 1 , E; 1 



-bi 

B n 




(8) 



-Bt 







Let us denote as z" . . . z^ m the 2m eigenvalues of T(e). 
The relations imply that if z™ is an eigenvalue of T(e), 
then (z~ n )* is an eigenvalue of T(e*). In this study we 
are concerned with the exponents 



£ a (e) = log|z a | 



(9) 



In general they still depend on n. Since | detT(e)| = 1, it 
is always ^ Q £ a (e) = 0. For real e the exponents of T(e) 
come in pairs ±£ a . 



C. Duality, Jensen, and spectrum of exponents 

Since the extremal components u\ and u n of the eigen- 
vector H(z n )u = eu enter in the eigenvalue equation @ 
of T(e), it follows that the characteristic polynomials of 
the two matrices are linked by a 

Duality relation, e is an eigenvalue of H(z n ) iff z n 
is eigenvalue o/T(e): 



det[e/ nm - H(z n )] (-l) r 



det[Bi •••£„] 



•det[T(e) - z n h 



(10) 



A proof of duality that holds also for non-Hcrmitian ma- 
trices, with blocks B\ being replaced by blocks Ck, is 
found in rcf. 4 . 

The spectrum of exponents can be obtained from the 
spectrum of H(z n ) through the following identity for an- 
alytic functions, which is a particular case of a theorem 
by Poisson and Jensen for subharmonic functions 10 : 

Jensen's identity. Let f be an analytic function in 
the open disk of radius R, where it has zeros z% , . . . , z^ 
that are ordered according to increasing modulus. Then, 
if < \z\\ and for r such that \zg\ < r < |z£+i| we have: 



/ W ^log|/(re^)|=log 
Jo 2vr 



^|/(0)| 

\Z! ■■■ Zt\ 



(11) 



Proposition. For real £ and complex e it is 

1 1 - 
- X)K-€a(e)]-e = ]Tlog|detB fe 



^ difi 1 

, , 2tt nm 



nm 



k=l 



log | det[e/„ m - 



)]| (12) 



Proof: Jensen's identity is applied to the polynomial 
f{z) = det[T(e) - z n I 2m ], with |/(0)| = 1 and z = 
gC+ji/V™. The duality relation is then used to obtain the 
formula. ■ 

For £ = a formula for the sum of positive exponents 
follows. It involves a real eigenvalue spectrum 

- E^( £ ) = -— E lo si dctBfe i 

m Z — t nm Z — / 

3 k 

dip 1 



«a>0 

2tt 



2n nm 



\og\det[eI nm -H(e^)}. (13) 



Equations (| 1 2[) and p3|) are exact and valid for a sin- 
gle, general transfer matrix. In the theory of disordered 
systems, a formula for the sum of exponents is known, 
where Jensen's angular average is replaced by the en- 
semble averag o n ' 12 i 13 . 

The left-hand side of eq. (fT!?]) is a non decreasing func- 
tion of £ (Fig[TJ) . For all £ > £max(c) (the maximum 
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FIG. 1: The behaviour of the right hand side of eq. (|12[l as a 
function of £, for m — 3. The constant value for £ < £i = £ m4n 
is + £2 + £3). At £ = £1,^2 and £3 = £a/ax the slope 
increases by 1/3. 



exponent of T(e)), the right side is always equal to £. 
For all positive £ < £ m m(e) (the smallest positive expo- 
nent), the right hand side is constant and equal to the 
average value of the exponents (|13p . For intermediate 
positive values of £ the function is piecewise linear, with 
discontinuities of order 1/m in the first derivative, at the 
values of the exponents. 



D. The matrix T f T 

Let us introduce the matrix Q(e) = T(e*) t T(e), with 
exponents r f a (e). For real e the matrix is real and posi- 
tive, and is preferred to T because of better large n be- 
haviour of the exponents. If B n B n = I m , the matrix Q 
is symplectic 



Q(e)E n Q(e) = £„ 



(14) 



and the exponents come in pairs ±7 Q . Hereafter, the 
matrix B n will be restricted to be unitary. Under this 
restriction, it is shown in Appendix A that (— l) n Q(e) is 
unitarily equivalent to the transfer matrix 9(e) 



'Ira " 




' I m " 


B n 


9(e) 


. B l . 



Q(e) = (-ir 

of the block tridiagonal matrix M(e) 
Ax — e Si 



(15) 



4-i 



B n -x 

A n — e I, 
I m e- 

B n 



A n Bj l _ 1 



B\ e-At 



Factors z ±2n are then introduced in the corners as in 
eq.© and, because of the factor (—1)™ in (fl3|) . the du- 
ality relation is as follows: 



dctM(e,z 2 ™) _ (-l) r 
IlJdetBfeP 



det[Q(e) - (iz) 2n I 2m ] 



./2nm, 



Since only the determinant matters, there is freedom to 
modify M to a form where e enters as a shift. Left 
and right multiplication by the block diagonal matrices 

\Inmi Imi Imi ~\~Imi • • •/ and {Inmi Imi Imi Imi ■ ■ ■} 

give detM = (~l) nm det[K((iz) 2n ) - el 2nm ], with 
K((iz) 2n ) = 

A 1 Bi 



Wr 



Bl 



B n -i 
T A 

B n -X 



-{izf 



••• Bl 

Bi A 1 



A true duality relation among eigenvalues is obtained 
(replace iz by z) 



det[K{z 2n ) - el 2nm ] _ (-1)' 



n fe |detB fc p 



y 2nrn 



det[Q(e) - z 2n I 2m ] (16) 



The Lyapunov spectrum is extracted by means of 
Jensen's formula: 

n-1 

£+- I>«(e)-e] = — E lo gi det ^i 



7a >« 



fc=l 



2tt 



dtp 1 
27r 2nm 



log|det[K(e 2 "«+^)-e/ 2nm ]|(17) 



Some properties of the matrix K are presented in Ap- 
pendix B. 



III. THE ANDERSON MODEL 

The discrete Anderson model describes a particle in 
a lattice, subject to a random potential. The poten- 
tial of a sample is specified by a set {vj} of random 
numbers chosen independently. Anderson^ considered a 
uniform density p(v) — 1/w in the interval [— w/2, w/2]. 
Lloy d 15 ' 16 studied the Cauchy distribution p(v) = £ (v 2 + 
<5 2 ) -1 , and evaluated the energy distribution exactly in 
any space dimension. Anderson's choice and the simple 
hypercubic geometry are here considered. More complex 
lattices can be studied by transfer matrix^. 

For a given configuration of potential, the eigenvalue 
equation is 



(18) 
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FIG. 2: Complex energy spectra for zero disorder, large n, 
£ = 1: in 2D (left, m = 5) and 3D (right, n x = n y = 4 i.e. 
m = 16). 



The sum is on the unit vectors along the 2D directions, 
e is the energy of the particle, the lattice has lengths 
rii, . . . ,tid- If the D axis is singled out, the sample is 
viewed as a number n = nu of sections each contain- 
ing m = m • • ■ riD-i sites. Accordingly, the Hamiltonian 
matrix is block tridiagonal 



H = 



Ai I n 



In 



J A 



(19) 



with Hcrmitian blocks Ai describing sections, and off di- 
agonal blocks describing hopping among sections. The 
associated transfer matrix is 



ne) = n 



el-n 



Ai 



-In 





(20) 



For large n the exponents of T(e) describe the inverse de- 
cay lengths of the eigenstates of Anderson's Hamiltonian. 
To study them, we introduce b.c. terms ±z™ in the cor- 
ner blocks of (|19p . and choose periodic b.c. in the other 
D — 1 directions, that appear in the diagonal blocks. 

Remark 0: For zero disorder the eigenvalues of H(z n ) 
are complex for any nonzero value of the parameter £ 
that measures non-Hermiticity (z = e^ +ttp ): 



Ree 



Im e 



2 cosh£ cos((/? 
2 sinh£ sin(<p - 



2tt 

n 
2ir 



i), 



£ = 1, . . . , n and r = 1, . . . , m. There are n eigenval- 
ues on each ellipse centered at e r = 2J2fS j 1 cos(2tt^-), 
1 < hi < rii- Therefore, the spectrum has support on 
m identical but shifted ellipses. In ID there is a single 
ellipse centered in the origin. In 2D there are m = n x 
distinct ones, while in 3D some of the m — n x n y ellipses 
may overlap because centers may be degenerate (Fig|5]). 
In Appendix C it is shown that the exponents of T and 
T^T coincide, for large n. 

Remark 1: For non-zero disorder the eigenvalues of 
the Hamiltonian matrix H(z n ) are all contained inside 
the union of ellipses 



(Ree - e ) 



(Ime) 



4 cosh £ 4 sinh £ 



< 1 



(21) 



where eo ranges in the interval [—2D + 2 — w/2, 2D — 2 + 
w/2]. Proof: If Hb(z)u = eu, and u is normalized, the 
inner product e = {u\H\ ) {z)u) in C nm is separated into 
real and imaginary parts: 

Ree - (u\Au) = 2|(n|S'it)| cosh£ cos{ip + 9) 
Ime = 2\(u\Su)\ sinh£ sin(<^ + 6) 

A is the block diagonal part of Hi,, S is the one-block shift 
matrix, and 9 =axg(u\Su) . The real number (m|^4u) = eo 
ranges in the spectrum of A. Schwartz's inequality gives 
the bounds. 

Remark 2: Since H is real, under complex conjugation 
it is T[e)* = T(e*). Then £ Q (e) = £„(e*). The symplectic 
property (|8|) with B n = I m implies that the exponents of 
T(e) come in pairs ±£ a for any e. 

Remark 3: Since the transposed matrix H{z n ) t coin- 
cides with H(z~ n ), then det[e/ nm — H(z n )] is a polyno- 
mial of degree m of the variable (z n + z~ n ). 

Remark 4- Since Hb{ze t2 ^^ n ) ~ Hb(z) (~ means sim- 
ilarity) and Hb(z)* = H^{z*), the following symmetry 
holds: H b (e^+^-^) ~ iJ 6 (e«+^)*. 



A. The Lyapunov spectrum 

The localization properties of Anderson's model are 
usually derived from the spectrum of positive Lyapunov 
exponents 71 < . . . < 7,„ of T(e) t T(e), with e real. Os- 
eledec's theorem^ guarantees that for large n it does not 
depend on the length n, and on the realization of dis- 
order. The most interesting exponent for physics is 71, 
that controls conductance. It is also the most difficult 
one to study numerically, because of the larger one s 19 ' 20 . 
Thorough investigations of the Lyapunov spectrum, its 
statistical properties and scaling, have been done in 2D 21 
and 3D—. The influence of b.c. was studied 2 ^ with the 
corner parameters z ±n of the present theory being both 
replaced by the same parameter t £ [0, 1]. It was found 
that the critical values of 71 and of the disorder parame- 
ter w c are t— dependent, while the critical exponent v is 
not. 

Analytic results for the Lyapunov spectrum are acces- 
sible in perturbation theory for the 2D strip 24 ' 25 i 26 , where 
n is large and m is finite. In such quasi-lD Anderson 
systems, the large n limit of the exponent spectrum of T 
coincides with the Lyapunov spectrum. For this reason 
here I concentrate on the spectral features of the matrix 
H{z n ) to which T is linked by duality. I'll then show in 
another section that the spectral identity (TT2")) allows to 
evaluate the smallest exponent £ m j n , which converges to 
71 for large n. 



B. ID Anderson model, Hatano and Nelson. 

For a chain of n sites the Hamiltonian is a tridiagonal 
n x n matrix and the transfer matrix is 2 x 2. Hatano 
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FIG. 3: ID Anderson model: Left: the complex eigenvalues 
form a closed loop, with two wings of real values (w = 7, 
n — 600, 5 = 1). Right: superposition of the eigenvalues of 
five matrices (same w and n) for £ from 0.5 (inner blob) to 1 
(outer). 



and Nelson^ suggested to study ID Anderson localization 
through the non-Hermitian extension of the model 



(22) 



with periodic b.c. As £ is increased from zero, the eigen- 
values do not distribute randomly in the complex plane 
but form a loop, Figj3] (left), whose analytic expression 
is known for Cauchy disorder- 2 ^.. The loop has two outer 
wings of real eigenvalues, that correspond to enough lo- 
calized eigenstates, and evolves to a more and more reg- 
ular shape, while the wings reduce, FigJ3] (right). The 
value £ c (e) up to which an eigenvalue e persists in the real 
axis, measures the inverse localization length of the phys- 
ical eigenvector (the Lyapunov exponent): £ c = 71(e). 
The latter is evaluated through Herbert, Jones and Thou- 
less' formula, 



71(e) = J deV(e')log|e-e'| 



(23) 



where p(e) is the disorder-averaged level density of the 
matrix ensemble in the limit of large n, £ = 0. The 
model has been studied by several author o 30 i 31 ' 32 i 33 ' 34 ; 
mathematical proofs were established by Goldsheid and 
Khoruzhenko 3 ^. Diagonalization of large non-Hermitian 
matrices is a delicate issue, as approximate eigenvalues 
may occur which are not close to true ones 3 -^ 

The Hatano-Nelson model is a case m = 1 of the theory 
presented in Section II. The duality relation flO]) simpli- 
fies greatly: 



det[e/„ - H(z n )} = trT(e) - (2 



and implies that 



trT(e) = det[e/ n - H(i)} = Pn {e) 



(24) 



(25) 



For a pure Bloch phase eq. (|24]) describes the energy 
bands of H(e mv ) as intersections of the polynomial y = 
Pn(e) with the strip y = 2cos(mp). As the non-Hermitian 
regime is entered, y = ±2cosh(n£), all eigenvalues of 
H{e n ^) are in the gaps, and approach pairwise for in- 
creasing £. A pair collides at a zero of p' n and becomes 
complex conjugate. This means that 2 cosh(nf ) equals 
the height \p n \ at an extremum of the polynomial. 



FIG. 4: ID Anderson model: the exponent £i(0) versus dis- 
order parameter w, for n = 600, averaged over 8 samples of 
disorder. 



The sum of the two eigenvalues e ± "^ 1+I¥ ' 1 ^ of the 
transfer matrix is 

p n (e) = cosh(n£i) cos(niy9i) + 2i sinh(n£i) sin(n^i) 

Elimination of the phase results in an exact equation for 
the exponent £i(e) 



(Rep n ) 2 



(Imp n ) 2 



4cosh 2 (n£i) 4sinh 2 (n£i' 



= 1 



(26) 



For large n it becomes e™^ 1 ^ = |p n (e)|, and gives a con- 
venient formula to compute the exponent (Fig[4]) 



log I det[e7 n - H(i) 



(27) 



In ID the single exponent is also given by the exact for- 
mula eq. (|13[) . For large n, £1 coincides with the Lyapunov 
exponent 71. 

Proposition: the eigenvalues of H(e n ^) distribute 
along the curve £i(e) = £. Real eigenvalues (wings) solve 
Pn(e) = 2cosh(n£). For large n£ the eigenvalues form the 
Lemniscat^Z- \p n (e) \ = e n ^. 



C. 2D Anderson model 

For a rectangular nxm lattice the Hamiltonian matrix 
(TT5]) has diagonal blocks 



Ai = 



Vi,l 1 

1 '•• 
1 



1 



1 Vi 



The eigenvalue spectrum of 2D non-Hermitian Anderson 
model is studied, and explained in the light of duality. 
Figj5] shows the eigenvalues of two matrices H(e n ^) with 
same £ and different m. They are distributed along a 
number of loops which is precisely given by m, the size 
of the blocks. 

By varying only the phase ip of z, the eigenvalues of 
the matrix Hb(z) move in the complex plane along arcs 
which retrace the loops. FigO shows that, as ip goes 
from to 27r/n, an eigenvalue moves along an arc that 



FIG. 5: 2D Anderson model: eigenvalues of a single matrix, 
with parameters w=7, n=100, £ = 1.5, ip = 0. Size of blocks: 
m = 3 and m = 10. The size m of blocks is the number of 
loops. 




FIG. 6: 2D Anderson model: motion of the eigenvalues in the 
complex plane for fixed disorder w=7 and parameters m—3, 
n=8, and £ = 1.5 and varying the phase (p. The 24 eigenvalues 
trace arcs: < ip < n /4 — S (left) and < <p < n/4 (right). 
The arcs join to form three loops. Loops are seen to contain 
different numbers of eigenvalues. 



FIG. 7: 2D Anderson model. Motion of 24 eigenvalues for 
< £ < 1.5, w=7, m — 3, n=8, (p = 0. The various wings 
terminate on m = 3 loops (at £ = 1.5 the phase is allowed to 
vary over 27T.) 




FIG. 8: 2D Anderson model: superposition of eigenvalue spec- 
tra for various realizations of disorder. For all: w=7, n — 100, 
£ = 1.5, <p — 0. Left: m=3, 20 realizations of disorder. Right: 
m=10, 5 realizations. The imaginary part of the eigenvalues 
is much less sensitive to disorder sampling than the real part. 



terminates where the arc of another eigenvalue starts. 
The union of such consecutive arcs makes a loop, and 
there are to closed loops. Differently from the w — 
case, loops may contain different numbers of eigenvalues 
of the matrix. The occurrence of loops is suggested by 
the duality equation: when a zero of dct[Hb(z) — el nm ] 
occurs, it is also a zero of det[T(e) — z n ], or 

2tt 

£a(e)=£, ¥>a(e) = ^ mod — , (a = 1, . . . , to) (28) 

n 

The loops are thus level curves of the exponents £ a , as 
functions of the complex variable e. In the limit of large 
n, the eigenvalues of one matrix fill m loops. 

In FigjT] only £ is varied: the eigenvalues trace lines 
that originate on the real axis (at £ = the matrix is 
Hermitian). For zero disorder the lines would be arcs 
of hyperbola. In the disordered case, for small £ the 
pattern of eigenvalues is complex, but evolves to regular 
loops. In the next section it is shown that the evaluation 
of exponents requires a "regular" regime. In Figj8] the 
parameters £ and ip are kept fixed, and the eigenvalues are 
computed for different realizations of disorder, with same 
strength w. They distribute along m loops, that appear 
shifted along the real axis for the different samples. 



D. The smallest exponent 

Let us assume that the exponents ±£ a of T(e) are iso- 
lated, < £ mm < £ 2 • • • < £m-i < £,max- By increasing 



£ from the value zero, the r.h.s. of eq. (fl2")) yields a con- 
stant value (the £ = value) until the value £ = £ TO i„(e) 
is reached. Then the function becomes linear with slope 
1/to until the value ^(e) is reached, where a new change 
of slope occurs. The change of slope can be used to 
identify the smallest exponent. FigJ5] illustrates this be- 
haviour. FigJTO] shows that for £ « £mm the eigenvalues 
of the matrix are all well in the complex plane. 

Kuwae and Taniguchi 3 - extended Hatano Nelson's ap- 
proach to 2D Anderson model, and evaluated numer- 
ically the average critical value £ c (tu) where the first 
pair of eigenvalues turns to complex. They conjectured 
that the inverse of the localization length coincides with 
this critical value, i.e. 71 = £ c . For example, for 
771 = n = 20, w — 7 they evaluate £ c « 0.1. Fig[S] 
shows that £, m in ~ 0.446, a measure of the inverse lo- 
calization length. For this £ the spectrum is evaluated 
and found in the complex plane. These facts adverse 
Kuwae's hypothesis that 71 signals the first escape to 
complex of the eigenvalues of Anderson's non- Hermitian 
Hamiltonian, and require further study. 

IV. CONCLUSIONS 

Based on a spectral duality relation for block tridiag- 
onal matrices and Jensen's identity, the distribution of 
exponents of a transfer matrix can be evaluated from 
the eigenvalue spectrum of the Hamiltonian with non- 
Hermitian boundary conditions. A preliminary numer- 
ical study of the complex energy spectra of Anderson 
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1.7164 
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FIG. 9: 2D Anderson model: evaluation of the r.h.s. of 
eg . (1 1 2 [I as a function of £ < 1, for e = 0, w = 7, aver- 
age on 40 angles. The change of slope marks £ m i n (0). Left: 
m = 3, n = 50, £ m in(0) « 0.87; the value 1.6692 is the av- 
erage + 6 + &)• Right: m = n = 20, £ mi „(0) « 0.447. 
The value 1.71627 is the average + • - • + £20). 




FIG. 10: 2D Anderson model: the eigenvalues for £ = 0.835, 
e = 0, iu = 7, ra = 3, n = 50. For £ near £ m i„ « 0.87 the 
inner loop is isolated. 



To obtain the structure (|29[) of a transfer matrix, it is 
necessary that B\B n = I m . The first and last factors 
containing B n are not consistent with (|29p . and only the 
intermediate product u\a\ ■ ■ ■ o\t\ is the transfer matrix 
of a tridiagonal block matrix, eq. (fl~5|) . 



Appendix B 

One can show the following properties of the matrix 
K(s). J is the matrix with 2n blocks I m along the di- 
agonal from lower left to upper right corners; S3 is the 
block diagonal matrix {/„,„, —Inm}- Then: 



JK{s)J = K{s*y 
S 3 K(s)S 3 = K(l/s*) 

The two imply that K(s) is similar to K(l/s) 



(30) 
(31) 



Appendix C 

For the Anderson model with no disorder (w = 0), 
real energy e and large n, the exponents of T and T^T 
coincide. Proof: It is 



Im 



(32) 



u " 




cl m A I m 


n 


'W 





u 




Im 







w 



non-Hermitian Hamiltonian matrices is made. The spec- 
tra have support on loops, that are explained as sections 
of the exponents at fixed height: £ a (e) = £■ This picture 
is complementary to a standard direct evaluation of the 
exponents £ a (e) at fixed energy e, by diagonalization of 
the transfer matrix. 

The spectral formula for exponents allows to evaluate 
the smallest one, and involves eigenvalues of the non- 
Hermitian Hamiltonian that are away from the real axis 
and distributed in loops that are nearly untertwined. 
This, and the higher exponents, are interesting subjects 
for further investigation. 

Appendix A 

We show the relationship of Q(e) = T(e*) t T(e) with a 
transfer matrix. It is convenient to factor T(e) as 



Im 

Accordingly: 

= (-1)" 



tn&n-l ■ • • <J\t\ 



Im 

Bt 



(29) 



Im 

B n 



Ui<t[ ■ ■■a] l _ 1 u n 



tk 



(BnBt)- 1 

I m 



dm ^4fe Im 
Ira 



Im 

Bl 



, u k = 

B^ 1 
B\ 



j4fe tlm Im 
Im. 



where A = UAU^ and A is the diagonal matrix of eigen- 
values {Ai, . . . A m }. The eigenvalues of T are m pairs 



zt n , where z k is a root of the equation z\ — (e— Xk)zk + 1 — 
0. 

The power n of the matrix can be computed by means 
of Cayley-Hamilton's formula. Because the blocks are 
diagonal, only powers zero and one of the matrix are 
needed: 





n 


' a " 




' /3 0" 




tlm A Ijyi 


Ira 




a 


+ 


(3 




Ira 



a and (3 are diagonal matrices with elements constructed 



with the roots Zk- 



Since e is real, cvfc 



and (3k are real. The matrix T^T is then constructed, and 
diagonalized. Its eigenvalues are pairs w^ 1 , with sum 



(*fc B + ** 2n ) 



zt + 1 



(33) 



If \zk\ > 1 then, for large n, \wk\ ~ \zk\ (the spec- 
trum of Lyapunov exponents of T^T and the spectrum 
of exponents of T coincide) ■ 
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